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Description 

[Title of the Invention] DIFFRACTIVE OPTICAL ELEMENT AND LASER MACHINING 
APPARATUS USING SAME 
[Technical Field] 
[0001] 

This invention relates to a diffractive optical element allowing a laser beam to branch 
into a plurality of beams and to shoot at any spots deviating from lattice points regularly 
assigned to an image plane and a laser machining apparatus using the same diffiactive optical 
element, particularly, a multipoint batch laser processing apparatus used for microprocessing of 
electronic parts, such as perforating, welding, soldering or surface treatment. 
[Background Art] 
[0002] 

[Patent Reference Document 1] Japanese Patent Publication No.2000-231012, "Method for 
designing a diffractive optical element" 

This is one of the former inventions made by tiie inventor of the present invention. The 
application proposed a diffractive optical elements (DOE: Fraunhofer type) having functions of 
dividing a laser beam. The mvention aimed at facilitating the designing of DOEs by assigning 
inherently narrow tolerance parameters witi. larger tolerances enlarged by the estimating 
designed DOEs by a merit function including influences of errors for the selected parameters. 

Enlargement of tiie tolerances makes tiie production of the DOEs easy. 

[0003] 

The DOE is a set of lengthwise and crosswise aligning identical unit patterns. A unit 

pattern is a set of lengflxwise and crosswise aligning pixels which are the smallest cell having a 

parameter. Since all the unit patterns are identical, the design of a DOE can be reduced to 

design of a unit pattern. The DOE with unit patten, diffracts a beam into a plurality of divided 
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beams just on lattice points assumed on an image plane. The Fast Fourier Transformation (FFT) 
^gorithm can calculate angular intensity distribution of diffracted beams separated by a DOE. 
[0004] 

[Patent Reference Document 2] Japanese Patent No.3,346,374, "Laser boring apparatus" 

This is another of the former inventions made by the inventor of the present invention. 
The patent No.3,346,374 proposed a laser boring apparatus for boring many holes with the same 
diameter at the same interval simultaneously on an object plane by dividing a laser beam into 
many identical separated beams and converging the separated beams by an fsinO lens in parallel 
on lattice points assumed on the object plane. Coupling of a DOE (Fraunhofer type or Fresnel 
type) and an fsinO lens raises position accuracy of the hole boring by the laser beams. The 
patent allows the apparatus to utilize higher order diffraction beams effectively without loss. 
The patented invention succeeds in boring many holes at a stroke and has various advantages. 
[0005] 

In the prior invention, the DOE was defined as an optical part which includes spatial 
repetitions of identical unit patterns T with a definite spatial period (interval) A and diffracts a 
laser beam into many divided beams. The DOE is made by designing a unit pattern T, aligning 
H identical unit patterns T in succession in a horizontal (x-) direction, and aligning G identical 
unit patterns T in a vertical (y-) direction. 
[0006] 

Then design of the DOE was equivalent to design of the unit pattern T . The DOE has 

a set of vertically-aligning R pixels and horizontally-aligning S pixels. A pixel is a minimum 

part whose thickness (heights) can take one of g=r (s: integer) different values. The full span 

of thickness difference of the pixels corresponds to a wavelength X of a laser beam. Pixels C are 

numbered by a horizontal number m and a vertical number n (0<m^H, O^n^G) as C{(m, n)}. 

Determination of a set of thicknesses {d„„} ofT different values for all the pixels C{(m, n)} was 
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design of a DOE. 
■ [O007] 

In practice, fte design of Ae DOE was far more aimplifiad. In fte prior DOE, all fte 
pixels could „„. take a arbiixary value of Sicknesses but were ruled by a r«ric,io„ condition. 
The r^swcaon is fte existence of a unit pattern T. The DOE had been designed by 
determining thicknes^s (d™) of a set of pixels ((m. n)} included in the unit pattern T and by 
ananging repeatedly GxH identical unit patterns in x- and y-directions. The design of the DOE 
with very high degree of fleedom could be replaced by the design of a unit pattern wifl, 1„„ 
degree of freedom. 
[0008J 

P x-direction aligning pixels and Q y-direction aligning pixels compose a unit pattern T. 
The unit pattern contains PQ pixels. The number of pixels as an object of design is not the total 
number GHPQ of total pixels but the number PQ of pixels within a unit pattern. 
[0009] 

A DOE has G x-aligning unit patterns and H y-aligning unit patterns. A unit pattern 
contains P x-aligning pixels and Q y-aligning pixels. A DOE has two-storied hierarchy. 
Parameters are defined here as follows. 
[0010] 

Pix€l C width (x-direction size) =a, length (y-direction size)=b. 
Unit pattern T PxQ pixels, width (x-size)=Pa, length (y-size)=Qb. 
DOE GxH unit patterns, width (x-size)^GPa, length (y-size)=HQb. 
DOE RxS pixels (R=PG, S=HQ) 

m, n coordinates of a pixel C in a unit pattern in the prior art 
dmn thickness of a C(m, n) pixel taking one of g=2^ values, 
tmn complex transmittance of a C(m, n) pixel. 
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m, n coordinates of a pixel C in the whole DOE in the present invention. , 

f 

P, q diffraction order parameters (integers) of the prior art 

a, P diffraction angles (not integers; continual) of the present invention. 
[0011] 

The set of notations is used throughout this description. A DOE includes RS (=PGHQ) 
pixels. A pixel takes one of the g=T (s: integer) different thicknesses. The total degree of 
freedom of a DOE is gRS(=2^PGHQ). But, actual degree of freedom of a prior art DOE is only 
gPQ, since repetitions of identical unit patterns compose a prior art DOE. The practical degree 
of freedom for designing is equal to gPQ which is a quotient of the total freedom degree gPGHQ 
divided by the unit pattern number GH. What is important is the number of pixels which should 
be detennined for designing a DOE. The pixel thickness parameter freedom g is less important 
here. Then, the degree of freedom for designing a DOE will signify the number of outstanding 
pixels (PGHQ or PQ) hereafter. 
[0012] 

For example, assuming that a unit pattern T includes 1000 pixels (PQ=1000) and a 
DOE includes 1000 unit patterns (GH=1000), the DOE includes 1000000 pixels 
(PQGH=1 000000). However, it is unnecessary to determine thicknesses of all the 1000000 
pixels. Instead, DOE designing should determine only the thicknesses of 1000 pixels. Free 
parameters are greatly reduced to one thousandth. 
[0013] 

Existence of repeating unit patterns alleviates the quantity of calculation to a great 
extent. The repetition of identical unit patterns gives a DOE a role of a diffraction grating which 
allows divided light to focus only at lattice points defined on an image plane. Extra regions 
except lattice points are all dark. The unit patterned DOE makes only beam spots only on lattice 

points which are distributed at predetermined spatial intervals (periods) on an image plane. The 
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lattice point diffraction enables the Fast Fourier Transformation (FFT) to calculate the intensity 
oYdiffracted beams at lattice points. Use of the FFT has exquisite advantages of fast calculation 
of a diffracted image and facile designing of the DOE. 
[0014] 

Furthermore, hole-boring processing, in many cases, bores many holes at lattice points 
alignmg lengthwise and crosswise at a common interval. For example, print circuit boards 
sometimes require that small holes are bored to be regularly arranged at a spatial period in two 
dimensions. 
[0015] 

Boring holes at lattice points fortunately corresponds to Fourier transfonnation of 
repetition of identical unit patterns. The FFT is a suitable calculation of the lattice point boring. 
Strictly speaking, since sine functions (sinO) of diffraction angles 9 are multiples of a unit value, 
an fsine lens instead of an ordinary ftanO lens allows the DOE to make lattice point diffraction. 

The lattice point diffraction can be automatically and rapidly calculated by the Fast Fourier 

Transformation (FFT). 

[Disclosure of Invention] 

[Problems to be Solved] 

[0016] 

Such a diffractive optical element (DOE) consisting of repetition of the unit patterns has 
drawbacks. 
[0017] 

Drawback (1): All that the DOE having equidistantly distributing unit patterns can do 

is to make divided beams only m regularly-distanced discrete directions leading to lattice points 

on an image plane. It is possible to shorten the pitch of a lattice by enlarging the size of unit 

patterns or by shortening a processing light wavelength X. The light wavelength is uniquely 
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determined by the kind of laser. Shortening of the wavelength requires replacement of a laser 
apparatus. 

[0018] 

However, we have not many kinds of high-power laser. High-power, available lasers 
suitable for laser processing are CO, lasers and YAG lasers at present. The processing hght 
wavelength should be otherwise determined by the material and the kind of processing (boring, 
annealing, welding). Thus, the wavelength should be intact. The wavelength is not a 
changeable parameter. Thus, reduction of the pitch of boring holes requires enlargement of unit 
pattems. 

[0019] 

A denotes a size (A=aP or bQ) of a unit pattern T. ^ is a wavelength of a laser. 
Integer j denotes the order number of diffracted beams. Integer m, means a dif&action order of 
the j-th beam, a, is a diffraction angle of the j-th diffracted beam. Bragg's diffraction formula 
holds, 

[0020] 

Asinaj=mjX,. (1) 

[0021] 

The diffraction order mj is a positive or negative integer (0, ±1, ±2, . . .). The left side of 
Eq.(l) is a difference between distances from two corresponding points on two neighboring 
pattems to an oblique plane slanting at an angle aj. The right side signifies a multiple of the 
wavelength. Eq.(l) means that the distance difference should be a multiple of the wavelength 
for diffracting light in the direction of Oj. 
[0022] 

An fsine lens enables the unit pattern-having DOE to converge diffracted beams at one- 
or two-dimensionally equidistantly distributed lattice points on an object plane (image plane) of. 



for example, a print circuit board. The lattice points on the object plane are distanced from the 
■ o-rigin by {fsina,} Q-O. «, ...). The common pitch p of the equidistant lattice points on the 
object is given by, 
[0023] 

p=f)JA. (2) 

[0024] 

Here, A is the size of a unit pattern; A=aP or bQ. Microprocessing requires further 
reduction of the pitch p. Eq.(2) signifies that an increase of the size A of the unit pattern will 
reduce the pitch p. 
[0025] 

The increment of the unit pattern size A, however, is restricted by some factors. The 
unit pattern size camiot exceed the size of the DOE. The lengthwise and crosswise repetition of 
identical unit patterns is indispensable for occurrence of diffraction. The quotient of the DOE 
size divided by the unit pattern size A should be a large integer. The total size of the DOE is 
i^stricted by the size (diameter) of an input laser beam. The laser beam diameter restricts the 
unit pattern size A via the DOE size. Thus, the maximum of A is determined by the input laser 
beam size. Thus, the pitch p should have a minimum determined by the maximum A. The 
existence of the large minimum p is an obstacle of further progress of microprocessing by the 
unit-pattern type DOE. 
[0026] 

Drawback (2): The conventional laser processing apparatus utilizing prior unit- 
patterned DOES makes dif&acted beam spots only on regular lattice points with a definite lattice 
constant on an object (image) plane but cannot produce diffraction spots on arbitrary non-lattice 
points on the object plane. El^tronic parts sometimes require the laser processing apparatus to 
make beam spots on sophisticated, non-lattice points on an object. The fsinO lens converges the 
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aj-diffracted beam on the mj-th point on an object, which is one of regularly aligning lattice 
points, 

[0027] 

fsinoj = mjXf/A. (3) 

[0028] 

Here, mj is an integer, f is a focal length, A is a unit pattern size, X is a wavelength and 
thus mjkf/A is the j-th lattice point defined on the object. 
[0029] 

The conventional unit pattern type diffractive optical element (DOE) is described by 
referring to Fig. 1. A DOE has GxR unit patterns T. A unit pattern T has PxQ pixels C. A 
pixel C has a length a and a width b (axb). 
[0030] 

Namely, PQ pixels build a unit pattern T. GH unit patterns construct a DOE. The 
DOE includes PQGH pixels. GH unit patterns have all identical distributions of pixels. Thus, 
degree of freedom of the DOE is only PQ, which is equal to the number of pixels contained in an 
T. 

[0031] 

When a plane wave light shoots the DOE, repetition of the unit patterns acts as a three- 
dimensional diffraction grating and makes a plurality of diffraction beams at the diffraction 
angles obeying Bragg's diffraction formula, i.e., Eq.(l) or El.(3). A horizontal diffraction order 
is denoted by p (integer) and a vertical diffraction order (y-direction) is denoted by q (integer). 
The total number of directions of diffraction beams is P in the horizontal (x-) direction and Q in 
the vertical (y-) direction. The right-hand figure in Fig.l is an object (image) plane T having a 
PxQ lattice structure. The lattice points coincide with crossing points of PxQ diffraction beams 

on the image plane. Any one mj of the PxQ lattice points defined on the object (image) plane T 
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is denoted by a set of integers (p, q) as order numbers. The origin is defined at the center point 
^ the object (image). The order numbers p and q satisfy inequalities of -P/2^p^P/2 and - 
Q/2^q^Q/2. 
[0032] 

The unit pattern included in the DOE takes another set of parameters (m, n) (m, n: 
integers) for numbering the PQ pixels. The x-direction (horizontal) order parameter is "m". 
The y-direction (vertical) order parameter is "n". The origin is determined at the center of a unit 
pattern. The unit order parameters m and n satisfy inequalities of .P/2^m^P/2 and - 
Q/2^n^Q/2. The complex amplitude transmittance of the (m. n)-th pixel in a unit pattern is 
designated by t„,„. In an ideal DOE, the absolute value | t.„ | of the complex transmittance t.„ 
isl. Thus, the transmittance can be written as t„„ = expOcp. J, where c^^^ Thephase 
cp.n is important, which is given by a phase difference between a phase change of vacuum 
passing and a phase change of pixel passing. 
[0033] 

The refractive index ofthe DOE is denoted by "n". The phase difference A<p is related 
to the thickness difference Ad„,„ ofthe m, n-th pixel by A<p.„ = 2.Ad„„(„.l)/X= kAd.„(n-l). 
Thickness steps (differences) Ad^ take g = 2^ (s; integer) different values. Phase differences 
Acp„„ and complex amplitude transmittance t.„ = exp(-jA<p.„)= exp(-j2;.Ad.„(n-l)) take g = 2'' 
different values. 
[0034] 

Although pixels are classified by black and white in Fig. 1, 2^ different steps (16, 32, 64, 
128. 256, 5 12. . .) are given to pixels in a practical case. When the DOE is irradiated with a laser 
beam, the DOE produces several diffraction beams in different directions. The diffracted beams 
converge at selected lattice points (p, q) on the object (image) plane T. Lattice points are 
defined as crossing points of horizontal lines and vertical lines. In any case, no beam is 
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diffracted to none-lattice point regions, which are defined as extra regions other than the lattice 

»• ^ 

points. Lattice pomts have possibility of being chosen to be irradiated. All the (PQ) lattice 
points, however, are not always irradiated. Some of the PQ lattice points are chosen to be 
irradiated but other lattice points are left unirradiated. The choice depends upon the unit pattern. 
[0035] 

The wording "lattice points" will be used also hereafter in this specification. Lattice 
points mean lengthwise and crosswise equidistant aligning points on an image (object) plane. 
Lattice points also denote equiangular points in diffraction angle space before shooting the image 
plane (T). The beams which pass the lattice points in the diffraction angle space will shoot the 
corresponding lattice points on the object plane (T). Directions or positions of the lattice points 
are exactly given by Eq. 2 or Eq. 3. 
[0036] 

Amplitude Wpq of the beam diffracted in the (p, q) direction by the DOE should be 
clarified. A lattice point (p, q) on an object corresponds to the diffraction angle order (p, q) one 
by one. Fourier transformation of the complex transmittance t,„„ of the DOE gives the (p, q) 
direction diffraction amplitude Wp,. A square | Wp J ^ of the complex amplitude Wp, is 
diffraction power at the (p, q) point on the object (T). 
[0037] 



[0038] 

Eq.(4) for Wp, is a simple equation suitable for calculation. Eq.(4) will be explained 

later. The summations SS of m and n should be carried not in the whole DOE but in only a 

single unit pattern T . The calculation is not a sum of PQGH pixels but a sum of PQ pixels. 

Integers p and q mean lattice points. Although it may not clear from Eq.(4), the summations of 
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W„ of E,.(4) for non-integers p and q are always zero. Integers p and , can yield non-zero W„. 
• Naniely. di£frac«on occurs only in fte dir«=tions (p. q) for integers p and q. Diffl^tion auZ 
beam spots only on the lattice points on an image plane. No difflaction beams make spots on 
non-latUce points. H^en. W„ should be calculated only for integet. p and q or only on lattice 
points (integer (p, q)). 
[0039] 

Calculation of Wp, for non-integers p and q is unnecessary, because it is known that 
Wp,=0 for non-integers p. q. This does not mean that an actual calculation of Wp, for non- 
integers p and q is always zero but the repetition of unit-patterns assures Wp,=0 for non-integers 
P and q. Then, a p, q non-integer Wp, should be left out of consideration. 
[0040] 

The above sine function is defined as sinc(x)=sin(;rx)/7ix which takes a maximum 1 at 
x=0 and oscillatingly decreases as x increases or decreases from 0. The inclusion of the 
function signifies reduction of higher order diffraction beams due to regular repetition of 
patterns in a DOE. 
[0041] 

When values of a set of {t„,„} are given, the Fast Fourier Transformation (FFT) 
calculate double summation EE in a short time without difficulty. Since the parameters p and 
included in expQ are integer., the FFT can calculate Wp, (Eq.(4)). The parameters p and q 
diffraction order parameters which correspond to mj in the Bragg diffraction formula Eq.(l). 
Thus, p and q are integers. The FFT can be only applied to the functions which have discrete 
parameters with a defmite sampling interval. The Fourier transformed function by the FFT is 
also has discrete parameters with another definite sampling mterval. Eq.(4) has the discrete 
transmittance matrix {t.„>. Wp, has discrete parameters p and q. Discreteness of t.„, p and q 
enables the FFT to calculate Wp, at a high speed. 
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[0042] 

The number of the parameters (m. n) on the DOE is PQ. The number of the parameters 
(P, q) on the image plane is PQ. The number of the DOE parameters is equal to the number of 
the image parameters. Beam power intensity I(p, q) in a (p, q) direction is given by squaring the 
absolute value of the complex diffraction amplitude Wpq. 
[0043] 

I(P,q)= I Wpq M (5) 

[0044] 

The design of a DOE is equivalent to the determination of Wp,. Wp, should be 
determined by presuming a desired diffraction beam pattern on an image plane, assuming t.„ for 
pixels on a DOE for producing the desired diffraction pattern, calculating Wp, from the assumed 
trnn by the Fast Fourier Transformation, squaring Wp, into I(p, q), comparing I(p, q) with the 
desired diffraction pattern, and examining the similarity. A plurality of candidates {t„„}s are 
obtained in general. An optimum {t^} can be chosen by estimating the candidates {t.„}s by a 
merit function. 
[0045] 

U, and V, are defined as two-dimensional angular resolutions of the conventional DOEs 
having repetition of unit patterns. Resolution means a minimum of discenuble spacing. 
Angular resolution signifies a minimum discernible angle. The angle between lines directed 
from the DOE to neighboring image lattice points is the necessary angular resolution for the 
conventional DOEs. The angular resolution for the conventional DOEs is e/L, where L is the 
distance from the DOE to the image plane and e is a lattice constant (= unit length of a lattice). 
Bragg-s diffraction condition is denoted by sina= mVA (m: integer, a: diffiaction angle). One 
step of a variation of sina is A(sina) =X/A. One step difference should be equal to the angular 
resolution. Then. e/L = X/A. The necessary angular resolutions U. and V. are a quotient of 
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a wavelength divided by the unit pattern size A(=aP or bQ). 
■ [0046] 

U| = X/aP. (6) 
V, =A,/bQ. (7) 

[0047] 

A DOE is an assembly of GxH unit patterns. A unit pattern is an assembly of PxQ 
pixels (axb). The width of the DOE is aGP=aR. The length of the DOE is bHQ=bS. The 
angular resolutions of the conventional DOEs camiot decrease less than a quotient (^aR or VbS) 
of the wavelength divided by the DOE size. The angular resolutions (minimum discernible 
angles) are too large in prior DOEs. 
[0048] 

Ui >X/aR. (8) 
V, >X^hS. (9) 

[0049] 

The rough resolutions mean that the conventional DOEs can make diffraction beams 
only upon lattice points and camiot diffract beams on non-lattice points. Object patterns must 
consist of lattice point spots. The confmed scope of the object patterns degrades the utility of 
DOEs. 

[0050] 

Derivation of Eq.(4) is briefly described. Wave optics requires to calculate image 
amplitude Wp, on an image plane (X, Y) by summing up individual amplitudes t„„ multiplied by 
oscillating parts given by expQ. The oscillating term expQ includes a part such as 
[0051] 

k(xX+yY)/L. (10) 

[0052] 
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^ L is a distance between the DOE and the image plane. Small letters x and y are 
• coordinates on the DOE. Capital letters X and Y are coordinates on the image. A size of a 
pixel, which has one definite parameter t.„. is axb. The coordinates x and y on the DOE can be 
converted into discrete representation; 
[0053] 

x=ma, y=nb. (11) 

[0054] 

Similarly, an image plane has a lattice structure with a lattice constant e. The 
coordinates X and Y on the image can be converted into discrete representation; 
[0055] 

X=pe, Y=qe. (12) 

[0056] 

Substitution of Eqs.(l 1) and (12) into Eq.(lO) produces 

[0057] 

k(mape + nbqe)/L. (13) 

[0058] 

An aperture angle looking from the DOE to the lattice length on the image is e/L, since 
e is the lattice constant and L is the distance fi-om the DOE to the image plane. A is employed as 
a symbol of discrete differentiation. Bragg's diffraction condition is expressed by e/L = A(sin a) 
= A(mjA/A) =X/A. 

The wave number k is given by k = 27r/X. The unit patter size is A=Pa or Qb. Eq.(13) 
can be converted into, 
[0059] 

(27oa)(map + nbq)X/A= 27t(mp/P+nq/Q). (14) 

The right-hand side of Eq.(14) is the same as the value included in the bracket of the 
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exp(,..) of Eq.(4). The above is a simple proof why E<,.(4) can give fte transformaUon from toe 
■ DOE paramete,.(w.o axe image pattern <W„>. Since m and p are integers from -P/2 ,o +P/2 
and n and p are integer, fton, -Q/2 .o+Q/2, E,.(4) having such an exp(...) can be swiffly 
calculated by the Fast Fourier Transformation (FFT). 
[0060] 

[Means for Solving the Problems] 

The present invention calculates the intensity W(a, P) of beams diffracted on non-lattice 
points (a, P) by assigning arbitrary complex transmittance {t.„} to all pixels (m, n) in a non-unit 
pattern DOE and summing up the following equation (15) discretely with m and n. 
[0061] 

Hc.fihsinc[fync[fj^P„ e.p{-Mn.aa.„l,fl)] (,5, 
[0062] 

The FFT has been partially abandoned by the present invention. Eq.(15) seems to 
resemble the FFT-suitable Eq.(4) in the form. However, Eq.(15) is entirely different from the 
FFT-available Eq.(4). Summation by m and n is not restricted within a unit pattern but done for 
all the pixels. Quantity of calculation increases up to GH times as much as the in-unit-pattem 
summation of the conventional FFT. 
[0063] 

In addition, diffraction point parameters are not integers p and q but non-integer 
fractional points which are denoted by continual angles a and p. (a, P) is not a lattice point on 
an image but a non-lattice point, a and p are not integers which prevent us from relying upon 
the FFT and force us to calculate the Fourier Transformation Eq.(15) by actually summing up 
with m and n step by step for all the pixels. 
[0064] 
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Besides the increment of the calculation quantity, the present invention increases the 
■ degree of freedom of designing DOEs. The increase of design freedom enables the present 
invention to assign diffracted beams upon non-lattice points on an image and enhance angular 
resolutions. But. the number of the diffracted beams (a, P) is still K which is equal to the 
number of the diffracted beams in the conventional FFT-operative DOE. The number of the 
object points does not increase. Inapplicability of the FFT and increment of tenns included in 
Eq.(15) up to HG times force the present invention to consume long time for calculation of 
Eq.(15) and assigning initial conditions and transmittance matrix {t„,„}. 
[0065] 

It takes long time to calculate Eq.(15). Then, the present invention uses Eq.(15) for the 
calculation of the intensity W(a„ p.) diffracted beams (a,, p J but uses the FFT for the 
calculation of noise. Noise, which should be zero in an ideal limit, is far smaller than W(a„ p.). 
The use of the FFT would bring about no serious error for the noise calculation. The present 
invention thus uses both the exact, direct summation calculation (Eq.(15)) and the convenient 
FFT calculation (Eq.(4)). This invention has a hybrid character. 
[0066] 

Fig.2 shows the relation between a DOE and an image of the present invention. The 
diffractive optical element (DOE) is not a set of repeated common unit patterns but a set of RS 
free pixels aligning in R lines and S columns, which can take free parameters. A conventional 
DOE has a three-storied structure of pixel/unit pattern/DOE. The present invention, which lacks 
the middle unit pattern, has a two-storied structure of pixel/DOE. 
[0067] 

A pixel can take any one of g steps of height (thickness), where g = 2, 4, 8, 16, 32. 
2' (s=I. 2. 3, ... ). As explained before. t™,=exp(-j2<n-l)Ad„7A.). where the height difference is 
Ad„n= {0, 1, 2, 3, (g-l)} X {A/(n-l)g}. 
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TTe DOE difflac a ia^r bean. i„.„ , bean,. i„ k di.c«o„.. Tie i„.e„.i., of ^ 
■ d.«.ac,ed in an , p. , ^ ^^^^ ^ ^ ^^^^^^ 

-nstaauon of a.e DOE co.p.ex .an.™i.^ce ,.„>. since *e OOE ine.udes no uni, 
pa«e™. U,ere is no concep. of .n-th orde. difiH,o.on any n,ore. Angular dependence of fte 
bean, difflac^d by ^ DOE can be ob^ned by casing on. aciuai c^cuiaUon of *e Fourie, 

transformation of {t,„„}. 
[0068] 

«'(«./») = sinc(i^jsi„.[Mj|;2'- exp{-,l(„„„.„,^), 
[0069] 

Here, k is Ae wave number 2^. What differs from E,.(4) is U,a. orien.a«on- 
designaiing paran,e.ers a and p are no. always integers (p and .0 but fractional «ai numbers in 
genera and U,a. *e scopes of sumnra«„ns of n, and n are no. a uni, pa..er„ bu, d,e whole DOB 
[0070] 

The change is not only a simple replacement from integer parameters p and q to real 
number parameters a and p. The change also means allowance of non-lattice angles which do 
not correspond to lattice points. Summations are not restricted within a small unit pattern (PxQ) 
but include all the pixels (RxS=PxQxGxH) in the whole DOE. 
[0071] 

Intensity distribution I(a, p) is obtained by squaring the complex amplitude W(a. P). 

[0072] 

I(a, p) = I W(a, p) I ^ (17) 

[0073] 

■■K" deno.es ,l,e number of tt,e diffi^ced. dive^g beams. K is a definite integer 

Which has been predetermined by the puroose of the nnp u 

y c puipose oi the DOE. K beams are designated by k=l , 2. 3, 
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...;K. Two-dimensional angles for the k-th diffracted beam are denoted by a. and fe. 
• [0074] 

Unlike the design of the conventional DOEs, it is not necessary to determine the 
diverging beam angles a and p as multiples of predetennined angular resolutions U, and V. (X/ 
A) in the design of DOEs of the present invention. However, it is still possible for the present 
invention to defxne necessary angular resolutions and V, other than X/A. Sines of the 
diffraction angles (a, P) of the diffracted beams can be determined as multiples of small common 
positive values U. and V. which should be properly chosen. Namely, sines of the angles can 
satisfy the following equations for some integers mk and nk, 
[0075] 

sinak = mkU2, (18) 
sinpk^n^Vg. (19) 

[0076] 

Here, k integers {mj should have no common divisor and k integers {n^} have no 
common divisor. Namely, (m J are prime to each other and K) are prime to each other. The 
greatest common divisor (measure) for groups {sinaj of sines of diffraction angles {a,> is 
defined as U,. The greatest common divisor (measure) for {sinpj should be defined as V, in 
the present invention. The multiples m^ and n^ are integers. Tl.us, U, and V, can be called 
angular resolutions for the DOEs of the present invention. 
[0077] 

Diffraction angles of the sine groups {sina^} and {sinPk} are not integers. The greatest 
common measures (GCMs) of (sina^} or (sinp^) may sound .trange. But, since {sina^} and 
{sinPk} are real number, the greatest common measure can be defined for the group numbers of 
{sinttk} or {sinpk}. 
[0078] 
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Here, a "common measure (divisor)" is defined as a non-integer number which can 
■ divide aU {sino,} wift integer quotients. The greatest common measu« is the greatest one 
among the common measures. The angle corresponding to Are greatest common measure is 
equal to the necessary angular resolution. 
[0079] 

Unless the angles and were small enough, inequalities of | sina,-sin a,- | <V2 
and I sinp,-sinp,. | <V, would hold. The required angular resolutions are given by Eqs.(18) 
and (19). 

[0080] 

Since a DOE consists of RxS pixels whose size is axb, the size of the DOE is aRxbS. 
All the pixels (m, n) have the freedom of taking arbitrary transmittance {t„,„}. The freedom of 
choice of {t„„} enables the angular resolutions Ua and V2 to satisfy inequalities U2<A/aR and V2 
<VhS. 
[0081] 

Here, the newly-defmed angular resolutions U2 and V2 are entirely different from the 
concept of the prior angular resolutions U,=AyaP and V,=X/bQ defined in the conventional design. 
The newly-introduced U2 and V2 are far smaller than the conventional U, and V,. 
[0082] 

One of the features of die present invention is that the DOE has no unit pattern T and 
all pixels are free parameters. All pixels can take an arbitrary, free value for the degree of 
heights. Since all the pixels are free parameters, the divergent beams can take free, arbitrary 
directions without the restriction of the Bragg condition. It is convenient to design the beam 
directions in free, arbitrary directions. However, the freedom of pixels deprives. DOEs of the 
regularity of unit patterns and the concept of diffraction orders. 
[0083] 
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The present invention, which lacks the concept of diffraction orders, should be still able 
' iS discern the minimum of differences of diffraction angles of diverging beams. The minimum 
angular resolutions are defined to be U2=AyaR and V2=X/bS. 
[0084] 

The newly-defined angular resolutions U2 and are still smaller than the conventional 
angular resolutions U. and V. (=VA) which are determined by the size of the unit pattern T. 
Increases of the number of parameters and the quantity of calculations can heighten the 
resolutions in the present invention. 
[0085] 

The present invention can be defined in different ways on different terms. This 
invention can be first defined as a beam diffraction apparatus which is founded upon high 
degrees of freedom in which all the RS(=PQGH) pixels can take free values of transmittance t.„ 
(m and n are line and column numbers of a pixel). The same matter can be expressed by 
inequalities of U2<XyaR and V^OJhS. The second definition based upon the inequalities may 
not be easily understood for the feature of the present invention. However, the second definition 
gives a criterion which is suitable for mathematical examination. 
[0086] 

The same matter can be otherwise represented by the request for extra calculations of 
diffraction intensity on non-lattice points other than the lattice points (p, q). However, the third 
definition of the present invention should be preceded by the definitions of diffraction images 
and lattice points of diffraction unages. 
[0087] 

Otherwise, the present invention can be negatively defined as an apparatus based upon a 
DOE which excludes repetition of similar unit patterns. 
[0088] 
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The present invention lacks unit patterns and pt^vents tlie Fast Fourier Transformation 
■ (FFT) on a unit pattern fh,m calculating an entire itnage pattern. Instead of flte patterned FFT, 
the present invention must calculate the complex amplitude W(a. P) at arbitrary, non-lattice, 
diverging angles a and 0 by summing up all Fourier transfom,ation contributions from all the 
individual pixels, "the full calculation of Fourier transformation takes long time in double 

meanings. The existence of die unit patterns restricts the angular resolution to X/A. The present 

mvention enjoys far smaller angular resolutions than VA. 

[0089] 

, The conventional unit patterned DOE makes difflBction spots of unit patterns only at 
lengthwise, crosswise, regularly aligned lattice points. No diffraction spots are produced on 
extra non-lattice regions. 
[0090] 

The present invention can obtain an arbitrary diffraction image having non-lattice 
diffraction points, since the present invention designs and makes a non-unit-pattem DOE. The 
number K of diffracted beams is definite and discrete. But the diffraction spots on an image 
plane are not always on lattice points. 
[0091] 

The image spots and the number of diffracted beams and the angular resolutions and 
V2 are contingent upon an actual pixel transmittance matrix {t.„}. The above multiple 
definitions of the present invention are not insignificant despite the contingency on {t„„}. 
[0092] 

It takes much time to calculate Fourier transfonnation for all pixels in the teaching of 
the present invention. An adequate selection of (t.„} enables this invention to bring diffraction 
beams at non-lattice points and obtain small resolutions and V. satisfying inequalities U,< 
X/aRand V2<A/bS. 
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[0093] 

•n.e conveaBonai uni.-panemed DOE canno, bring diffiaotion beams a. no„-laHice 
points and camio. reduce angular resoluUons U, and V, less man VaR and V.<VbS for any 
selection of {tmn}. 
[0094] 

It is not easy to understand the above differences between the conventional unit- 
patterned DOE and the present invention DOE. 
[0095] 

Exclusion of unit patterns and overall free choice of pixels allow the present invention 
to produce diverging beams aligning in arbitrary spot numbers and distributions. 
[Effects of the Invention] 
[0096] 

•nie present invention can irradiate a plurality of arbitrary points simultaneously with 
laser beams with high accuracy. This invention succeeds in free, precise laser simultaneous 
multipoint processing which is required for electronics products. Another advantages of the 
present invention are acceleration and cost-reduction of the laser processing. 
[0097] 

The conventional DOE, which has repetition of common unit patterns, makes a plurality 
of beams diffracted into predetermined, discrete lattice points aligning two-dimensionally at 
constant intervals (lattice constants). The function of the DOE is essentially Fourier 
transformation for making divided multibeams shooting upon the known lattice points regularly 
aligning on an object. Thus, it was enough for the conventional DOE to calculate the diffiaction 
beam intensity only on the predetermined object lattice points. 
[0098] 

The conventional DOE enables the Fast Fourier Transfonnation (FFT) to calculate 
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diffiactio. beam. i„ a short toe. However, U,e unit-patterned conventional DOE malces 
• diffraction beam, jus, on the predetermined lattice points a, intervals (lattice constants) which 
are equal to X/A. where X is a wavelength of a laser and A is the size of a unit pattern. Tk. 
lattice constant is pretty large, which is one drawback. A wide lattice constant is toconveniem to 
depict fine difSacted patterns. Reduction of the lattice constant requires enlargement of the unit 
pattern. But. the enlargement is not always possible, since the unit pattern size should be a 
divisor of a DOE size. Another dtawback is that it is impossible to produce non-lattice. free 
patterns includmg non-lattice pomts which deviate from lattice points. 
[0099] 

Instead of building up repetition of common unit patterns, this invention gives all pixels 
fbll fieedom of flrickness (or height). Thus, the present invention increases the degree of 
freedom up to GH times as much as the conventional DOE. where GH is the number of unit 
patterns moluded in the DOE. TTre points shot by diffraction beams are not always lattice points. 
[0100] 

The concept of lattice points vanishes in the present invention, which dispenses with 
unit patterns. The spot distribution changes from discrete, digital distribution to analog 
distribution. 

[0101] 

The number of points 6n which the diffraction intensities should be calculated increases 
to a great extent. The Fast Fourier Transformation is inapplicable. TTie difiraction intensity at 
necessary points should be calculated by the whole Fourier transformation. These are weak 
points of the present invention. 
[0102] 

Despite the drawbacks, the present invention has strong points of producing arbitrary 
patterns consisting of non-regular, non-lattice point distribution on an object, enhancing 
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precision by fining resolutions, widening the scope of applicability of DOEs, and heightening 

utility of DOEs. 

[0103] 

This invention enables the diffractive optical element to irradiate a plurality of arbitrary 
points simultaneously with high accuracy with divided laser beams. This invention satisfies a 
requirement of high-precision laser light processing at arbitrary points, which is suitable for 
making electronic devices. This invention is effective for acceleration and cost-alleviation of 
laser processing. 

[Best Mode for Carrying Out the Invention] 
[0104] 

The above calculation for arbitrary angle diffracted beams enables the present invention 
to optimize the design of DOE patterns. There are several optimizing methods. Here, a Direct 
Binary Search (DBS) method is adopted as an optimizing method like our prior designing of 
Japanese Patent Publication No.2000-231012 of the same inventor. 
[0105] 

The DBS method employs a merit function for estimating the characters of DOEs. The 
DOE patterns are optimized by giving trial (t„,„} to pixels, calculating W(a, P) for the {t„„} by 
Eq.(16), estimating the result by the merit function, changing {t.„} for the direction of lowering 
the merit function, and minimizing the value of the merit function. The set of {t„„} which 
minimizes the merit function is one of the candidate transmittance matrixes {t,„}. 
[0106] 

When a set of {t„,„} is given, W(a. P) can be calculated by Eq.(16). Beam intensity 
fmictions I(a, P) is given by I(a. P)= | W(a. P) | \ For example, the following merit function E 
allows us to estimate the solutions, 
[0107] 
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, f="'^('7o-En«».A)J + »'.2:(/(a.,A)-77 (20) 
[0108] 

In the expression, is an object intensity value of the diffraction beams, T is an average 
of I(a. P), and W. and are weight constants. An adequate set of {t„,„} is sought by 
minimizing the merit function E. In minimizing E. the first term of the right hand of Eq.(20) has 
a tendency of raising diffraction efficiency and the second term has an effect of enhancing 
uniformity of diffiracted beams. Minimizing enables the merit function to design a suitable DOE 
which diffracts a beam into desired directions with high diffraction efficiency and uniformity. 
[0109] 

The merit fianction of Eq.(20) neglects noise and reduction of diffracted beams by noise. 
The use of Eq.(20) would allow strong noise to occur in undesired directions. Besides the main, 
processing beams, hoise should be taken into account. However, calculation of noise is more 
difficult than the calculation of diffraction beams. Of course, noise can be calculated by W(a, P) 
of Eq.(16). But, noise happens in unexpected directions. Unlike the programmed diffi-action 
beams, the directions (a. P) of noise occurrence cannot be prophesied. To investigate noise 
occurrence in all dkections, all W(a, P) should be calculated except diffraction directions (a., p,). 
It is nearly impossible to calculate all the W(a, P) for continual (a, P) except diffraction (a^, p,). 
[0110] 

On the contrary, the conventional calculation based upon the FFT can calculate PxQ 
spectra including noise by a single cycle of calculation. 
[0111] 

Then, this invention calculates noise distribution by the Fast Fourier Transformation 
(FFT) which can obtain only Wp, at lattice points (p. q) (p. q=integers). The FFT cannot 
calculate noise at non-lattice points. But, it matters little. Noise occuixence is continual. Non- 
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lattice point noise can be approximated by the nearest neighboring lattice point noise. The 
• present invention has a hybrid character. Diffraction beams are calculated by the entire Fourier 
transformation calculation in the whole DOE, but noise calculation is carried out only at lattice 
points by the FFT as expressed by Eq.(4). A discrete angular spectrum which is obtained by the 
Fast Fourier Transformation (FFT) for the RxS pixels (all pixels in DOE) is denoted by I.. I. 
still has diffraction beam intensity as well as noise. I„ takes large values near the k-th 
diffracting beams (a,, p.). Then, IV. is defined as a value by subtracting diffraction beam 
intensity from the whole spectrum I„. Then, IV, includes only noise spectrum. For example, 
ir^stead of the merit function Eq.(20), noise can be taken into consideration by using the 
following new merit function, 
[0112] 

^ = ^{'7o-i:^k.^j] +^2i:fe./^J-7j+fF3Max(/'„y (21) 
[0113] 

The third term of the right side reveals a function of reducing noise (W3 is a weight 
constant). Noise calculation is prosecuted by the FFT. Thus, the noise angular resolutions are 
restricted by U3=VaR and V3=X/bS. 
[0114] 

Since the present invention essentially designs a DOE without using the FFT, there is a 
possibility that the noise may be included at narrow, localized regions below the large 
resolutions U3 and V3. The above FFT cannot detect such localized, hidden noise. The inventor 
further contrives two noise detection methods. 
[0115] 

[Noise detection method (1)] 

The calculation region is enlarged to R'xS' pixels (R'>R, S'>S) by adding virtual, 

26 



vacant pixels of U=0 to neighboring regions around the RxS pixels. Enlargement of the 
ciilculation range reduces resolutions down to W=yaR' and V3'=AybS'. Of course, U3'<U3 and 
V3'< V3. Since angular resolutions are refined for noise, noise detection becomes facile. 
[0116] 

[Noise detection method (2)] 

Eq.(16) does not include the information of incident beam intensity distribution. The 
present invention is also applicable to another design taking account of the incident beam 
intensity distribution. 
[0117] 



W{a , fi) = sinc[^]sinc[Mj2:Xa„„r„„ exp{- jk{n,aa nb/j)} 



(22) 



[0118] 

Here, a„,„ is complex amplitude of the incident beam. The beam is divided into discrete 
small rectangle cells which have the same size as the pixel of the DOE. Then, the divided beam 
cells are numbered by m and n which are the same as the pixels. The (m, n) cell of the incident 
beam aligns along the same axial line which pierces the (m, n) pixel of the DOE. 
[0119] 

For example, if the incident beam is a Gaussian beam with a l/e^ diameter of aR/2, an 
angular diameter of the Gaussian beam in an angular spectrum is SXXnaR = %U^Xn, which is 
over twice the afore-mentioned angular resolution U3. A noise spot has a similar diameter to the 
beam spot. This makes it easy to detect noise by the FFT irrespective of the restriction on the 
angular resolution. 
[0120] 

Eq. (16) aims at calculating an angular spectrum on an image, that is, a far-field pattern 
(Fraunhofer diffraction image). The DOE designed by the teaching expressed by Eq.(16) 
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Fourier type (=Fraunhofer type) DOE having an indefinite far-distanced focus (f = oo). 
■ [0121] 

A Fresnel type DOE having a definite focal distance f can be obtained by superposing 
the Fourier type DOE upon a pattern made by a converging lens with a definite focus f Before 
the superposition, the pixel size of the Fourier type DOE should be equalized to the pixel size of 
the lens pattern. 
[0122] 

A pattern of a shorter focal length lens requires smaller pixel size. If the focus of the 
lens is too short, it is preferable to divide a pixel of the Fourier type DOE from a x b into a' x b', 
where a'=a/C., b'=b/C2, (C., C^; integers). The converging lens should be designed on the 
reduced sizes a' and b'. Then, a Fresnel type DOE is made by superposing the Fourier DOE of a 
unit size of a' x b' on the lens of a unit size of a' x b'. 
[0123] 

A Fresnel type DOE can reduce the influence deriving from microprocessing errors (e.g., 
step errors), in particular, reduce the fluctuation of the 0-th order beam (a=p=0), which had been 
described in our Japanese Patent No.3,346,374. 

A Fresnel type DOE can be designed by multiplying a lens pattern to the Fourier type 
DOE of the present invention. Thus, the Fresnel type DOE's W(a, P) can be obtained by 
summation of an equation similar to Eq.(16) with m and n with regard to all pixels. 
[Embodiment 1] 
[012^1] 

[EMBODIMENT 1 (Fraunhofer type (= Fourier type) DOE, small diverging angles)] 
[A. Properties of the Fraunhofer DOE] 
Wavelength: 10.6 jxm 
Mode: TEMqo 
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Beam diameter(l/e^): 10mm 

Step number of heights: 16 

Pixel size: a = 20 }im, b = 20 nm 

Pixel number in a DOE: R = 1000 pixels, Q = 1000 pixels 

DOE size: aR = 20 mm, bS = 20 mm 

Number of divided beams: 49 (=7x7) 

Diffracted beam pattern: Fig. 3, P, = 3.93701 mrad, P2 = 3.93701 mrad, 

6, = 0.15748 mrad, 82 = 0.07874 mrad (lad = 57.29 degrees, 

mrad= 10~^rad) 
[0125] 

Angular resolutions U and V are sought by the equations sina,=m,U and sinp,=n,V. In 
Embodiment 1. U = 6, = 0.15748 mrad and V = 6. = 0.07874 mrad. X/aS = M>S = 0.53 mrad. 
Embodiment 1 satisfies the prescribed conditions U<yaR and V<A/bS. 
[0126] 

[B. Result of designing of DOE] 
DOE of Embodiment 1 

Diffraction rate: 82.2% 

Intensity fluctuation : 1 .6% 

Largest noise intensity: 3.6% 

Fig.4 is an incident beam intensity pattern. Fig.5 is a DOE step distribution pattern. 
Fig.6 is an image spot distribution pattern. 
[0127] 

[C. Laser processing apparatus] 

A laser processing apparatus is manufactured by combining a CO, laser, tl,e above- 
mentioned Fourier type DOE and a flat/convex lens (f - ,27 mm. 50.8 mm diameter. ZnSe). 
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The spot diameter on an image plane is 171 |xm(p. 

•I 

[0128] 

An object arrangement of spots in angular space of Embodiment 1 is shown in Fig.3. 
The abscissa is a horizontal angle a. The ordinate is a vertical angle p. It should be noted that 
this pattern is not an image pattern projected onto the image plane. 7x7 spots align in a- and p- 
directions. The beam spots are not all upon lattice points but slightly deviate from the lattice 
points. 

[0129] 

Four (±2, ±2) points, four (0, ±2), (±2, 0) points and one (0, 0) point are on the lattice 
points. Other points deviate from the lattice points. Vertical and horizontal series of the spots 
depict wavy lines instead of straight lines. The off-lattice pattern is the object of Embodiment 1. 
The conventional unit-patterned DOE cannot produce such an off-lattice pattern. 
[0130] 

Lattice constants are P,=500 irni and P2=500 pixn. A lattice has a size of 500 pmixSOO 
Mm. hi Fig. 3, four (±3, ±3) points deviate by 5, and 82 in a counterclockwise direction from 
lattice points. The (+3, +3) point is off from the nearest lattice point by -6,= -20 pxn and 83= 10 
^m. Deviations are exaggerated in Fig. 3 for emphasizing that the invention DOE can make off- 
lattice diffraction beams. Another (-3. -3), (-3, +3), (+3,-3) points deviate by ±6,(=20 Mm) and 
±52(=10 iim). Since given deviations are 6, and 62, resolution angles U3 and V3 should be 
smaller than 61 and 62. 
[0131] 

Such off-lattice diffracted beams cannot be produced by the conventional DOE with 
repetition of common unit patterns. What produces the off-lattice pattern is only the DOE which 
the present invention proposes. 
[0132] 
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Fig. 4 is a spatial power distribution of the incident laser beam. The incident beam is a 
CO2 laser beam of a Gaussian power distribution. 
[0133] 

Fig. 5 shows the pattern of the DOE designed by the teaching of the present invention. 
The DOE seems to be repetition of similar amorphous blank islands. But, the blank islands are 
different in shape and size. The DOE is not repetition of identical unit patterns. 
[0134] 

Fig. 6 is a 7x7 diffraction beam spot angular distribution on an image produced by the 
DOE of Embodiment 1. The size of the image including margins is 39.75 mrad. A neighboring 
spot interval is nearly 5mrad. Points are not on lattice points but deviate from virtual lattice 
pomts. Every set of vertical seven points is not on a direct line but on a wavy line. All third 
line points and all seventh line points deviate to the right. All first line points and fifth line 
points deviated to the left. The 7x7 object pattern includes diffraction points deviating from 
virtual lattice points. Embodiment 1 allows the beam spots to deviate arbitrarily from the lattice 
points. The beam intensity is uniform. 
[Embodiment 2] 
[0135] 

[EMBODIMENT 2 (Fraunhofer type (= Fourier type) DOE. large diverging angles)] 

[A. Properties of the Fraunhofer DOE] 

Wavelength: 10.6 jam 

Mode: TEMqo 

Beam diameter(l/e^): 10 mm 

Step number of heights: 1 6 

Pixel size : a = 20 jxm, b = 20 ^m 

Pixel number: R = lOOO pixels, S=1000 pixels 
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DOE size: aR = 20 mm, bS = 20 mm 

Number of divided beams: 49 (=7x7) 

Diffracted beam pattern: Fig. 3, P,= 39.37008 mrad, 39.37008 mrad. 

8i= 0.15748 mrad, 62= 0.07874 mrad 

[0136] 

The object diffracted beam pattern is similar to the pattern shown in Fig. 3 having 
fluctuating (7x7) spots deviating from lattice points for clarifying difference between 
Embodiments 1 and 2. The patterns are similar figures but the sizes are different between 
Embodiments 1 and 2. Embodiment 2 chooses lattice constants P,= 39.37008 mrad and P^= 
39.37008 mrad for investigating a case of large diverging angles. Embodiment 2 has the lattice 
constants which are ten times as large as Embodiment 1. Deviation of spots from the lattice 
points, however, are the same as Embodiment 1. Namely, Embodiment 2 takes 5,= 0.15748 
mrad and 5^= 0.07874 mrad which are the same as Embodiment 2. Ratios of deviation/lattice 
are reduced. 

[0137] 

In Embodiment 2, necessary resolutions U and V in x- and y-directions are 
U = 5] =0.15748 mrad. 
V = 62 = 0.07874 mrad. 

[0138] 

In Embodiment 2. aR = 20 mm (20000 jrni) and bS = 20 mm (20000 jim). Ratios of 
wavelength/DOE are VaR=X/bS=0.53 mrad. These values of parameters satisfy inequalities, 
[0139] 

U(0. 1 5748 mrad) < VaR(0.53 mrad), 
V(0.07874 mrad)<A/bS(0.53 mrad). 

[0140] . 
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Embodiment 2 also satisfies the prescribed conditions. 

' [0141] 

[B. Result of designing of DOE] 
DOE of Embodiment 2 
Diffraction rate: 71.6% 
Intensity fluctuation: 1 .7% 
Largest noise intensity: 3.5% 

The pattern of the DOE is depicted by Fig.7. Thickness (height) distribution (16 steps) 
of the DOE is represented by a black-white monochromatic tone. Fig.7 shows a lattice-like 
pattern with narrow lattice constants. A fine structure of tones within lattices are obscure. The 
DOE (Fig. 7) of Embodiment 2 is entirely different from the former DOE (Fig. 5) of 
Embodiment 1. Since the branching beam angles are large, repetition of the DOE pattern is fine 
in Embodiment 2. 
[0142] 

No unit pattern should exist in the DOE in Fig. 7. There should be no repetition of unit 
patterns having the same shape and the same size. But. the DOE seems to have quasi-regular 
repetition of spots in two dimensions like cloth. It may be marvelous. But, the deviations 6. 
and 62 are quite small in comparison to the pitches Pj and P2 in Embodiment 1. Then, the DOE 
becomes similar to the conventional regular DOEs consisting of repetition of unit patterns. 
[0143] 

The solution of designing the DOE which produces spot disposition of Fig. 3 is not 
unique. Patterns of the DOE for realizing an object spot distribution have many solutions due to 
many degrees of freedom of the DOE. A pertinent solution has been chosen by a merit function. 
The merit function selects a suitable pattern for reasons of a high diffraction ratio, satisfaction of 
object values or low fluctuation of power of branched beams. 
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[0144] 

[C. Laser processing apparatus] 

A laser processing apparatus is manufactured by combining a CO^ laser, the above- 
mentioned Fourier type DOE and an fsinG lens (f = 127 mm). The spot diameter on an image 
plane is 171 ^mcp. This fsin9 lens is produced by the teaching of the previously described 
Japanese Patent No.3,346,374 invented by the same inventor as the present invention. 
[0145] 

[Embodiment 3] 

[EMBODIMENT 3 (Fresnel type DOE, small diverging angles)] 
[A. Properties of the Fresnel DOE] 
Wavelength: 10.6 ^m 
Mode: TEMoo 

Beam diameter (1/e^): 10 mm 

Step number of heights: 16 

Pixel size: a = 20 fim, b = 20 ^m 

Pixel number: R = 1000 pixels, S = 1000 pixels 

DOE size: aR= 20 mm, bS= 20 mm 

Number of divided beams: 49 (=7x7) 

[0146] 

Diffracted beam pattern: Fig. 3, P,= 3.93701 mrad. 3.93701 mrad, 

5,= 0.15748 mrad, 5^= 0.07874 mrad 

[0147] 

Fresnel type lens: focal length f = -500 mm 
[0148] 

What is different from Embodiment 1 (Fraunhofer type) is that the DOE has a function 
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as a lens. Since the DOE has the lens function, the DOE is called as Fresnel type. Embodiment 
■ 3 adds an -500 mm concave lens to the DOE. Of course, a convex lens can be added to the 
DOE instead of the concave lens. 
[01 49] 

Necessary resolutions U and V in x- and y-directions are 
U =6, =0.15748 mrad, 
V =52 = 0.07874 mrad. 

[0150] 

In Embodiment 3, X/aR = A/bS = 0.53 mrad like Embodiment 1. These values of 
parameters of Embodiment 3 also satisfy the inequalities, 
U<AyaR, 
V<AybS. 

[0151] 

[B. Result of designing of DOE] 
DOE of Embodiment 3 

Diffraction rate: 82.2% 

Intensity fluctuation: 1 .6% 

Largest noise intensity: 3 .6% 
[0152] 

[C. Laser processing apparatus] 

A laser processing apparatus is manufactured by combining a CO^ laser, the above- 
mentioned Fresnel type DOE and a flat/convex lens (f = 127 mm. 50.8 mmq,, ZnSe). The spot 
diameter on an image plane is 171 nm(p. 
[0153] 

Fig. 9 is a height (thickness) pattem of the DOE of Embodiment 3 designed by the 
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teaching of the present invention. Amorphous islands forming parts of concentric circles are 
• prevalent. The concentricity derives from the lens character which is inherent in Fresnel type 
DOES. The reason why the island sizes are large is that the lattices P, = 3.93701 mrad and = 
3.93701 mrad are small. 
[Embodiment 4] 
[0154] 

[EMBODIMENT 4 (Fresnel type DOB, large diverging angles)] 

[A. Properties of the Fresnel DOE] 

Wavelength: 10.6^ 

Mode: TEMqo 

Beam diameter(l/e^): 10 mm 

Step number of heights: 1 6 

Pixel size : a = 20 jxm, b = 20 ^m 

Pixel number: R = 1 000 pixels, S = 1 000 pixels 

DOE size: aR = 20 mm, bS = 20 mm 

Number of divided beams: 49 (=7x7) 

Diffracted beam pattern: Fig. 3, P,= 39.37008 mrad, P^= 39.37008 mrad, 

6i= 0.15748 mrad, 82= 0.07874 mrad 

[0155] 

Fresnel type lens: focal length f = -500 mm 
[0156] 

An object pattern is a fluctuating 7x7 dot pattern shown in Fig. 3. Embodiment 4 takes 
P.= 39.37008 mrad and P,= 39.37008 mrad for investigating the influence of a large diverging 
angle case. Embodiment 4 is similar to Embodiment 2 in the large diverging angles. 
Embodiment 4 differs from Embodiment 2 at the point of the DOE having a combined lens. 
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[0157] 

Embodiment 4 makes a DOE shown in Fig. 10. The DOE of Fig. 10 is different from 
Fig. 7 (Embodiment 2). Fig. 7 shows small dots aligning lengthwise and crosswise. Fig. 10 
shows concentricity of dot arrangement. The concentric character derives from the lens 
combined. 
[0158] 

Necessary resolutions in x- and y- directions are 
U= 6, = 0.15748 mrad, 
V= 62 = 0.07874 mrad. 

[0159] 

In Embodiment 4, aR = 20 mm (20000 ^m), bS = 20 mm (20000 jmi) and X = 10.6 [im. 
Then, A/aR=AybS=0.53 mrad. These values of parameters of Embodiment 3 also satisfy the 
inequalities, 
[0160] 

U(0. 15748 mrad)<AyaR(0.53 mrad), 
V(0.07874 mrad)<XybS(0.53 mrad). 

[0161] 

[B. Result of designing of DOE] 
DOE of Embodiment 4 

Diffraction rate: 71.6% 

Intensity fluctuation: 1.7% 

Largest noise intensity: 3.5% 
[0162] 

[C. Laser processing apparatus] 

A laser processing apparatus is manufactured by combining a CO2 laser, the Fresnel 
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type DOE and an fsinG lens (f = 127 mm). A spot diameter on an image plane is 171 ^m<p. The 

fsinG lens has been produced by the teaching of our Japanese Patent No.3,346,374. 

[0163] 

Fig. 10 is a height (thickness) pattern of the Fresnel type DOE of Embodiment 4 
designed by the teaching of the present invention. Plenty of small dots forming parts of 
concentric circles are prevalent. The concentricity derives from the lens character which is 
inherent in Fresnel type DOEs. The reason why the dot sizes are small is that the lattices Pi= 
39.37008 mrad and P2= 39.37008 mrad are large. It is reasonable that Fig. 10 seems to be a 
superposition of a concentric lens pattern upon the Fig. 7 DOE pattern. 
[Brief Description of the Drawings] 
[Fig. 1] 

Fig. 1 is an explanatory figure of a conventional dififractive optical element (DOE) and 
an image plane, the DOE consisting of repetition of identical unit patterns with PxQ pixels, 
having a degree of freedom of PxQ, making diffracted spots only on PxQ lattice points 
periodically two-dimensionally defined on the image plane, and allowing the Fast Fourier 
Transformation calculation (FFT) to calculate diffraction light intensity distribution Ipq only at 
the PxQ lattice points. 
[Fig. 2] 

Fig. 2 is an explanatory figure of an inventive diffractive optical element (DOE) and an 
image plane of the present invention, the diffractive optical element having not unit pattern but 
freely-determinable RxS pixels which require to integrate a diffraction integration equation for 
calculating diffraction light intensity distribution on the image plane instead of the FFT and the 
image plane, and having no concept of lattice points. 
[Fig. 3] 

Fig. 3 is an angular distribution of forty-nine beams divided by the DOEs of 

38 



Embodiments 1 to 4 of the present invention, wherein beam spots are not exactly on the lattice 
points and deviate from lattice points. Deviations are exaggerated for attracting attention. 
Deviations 61 and 62, which are the minima of necessary resolutions, are far smaller than lattice 
intervals (lattice constants). The conventional method was satisfied with a resolution as large as 
a lattice constant. On the contrary, the present invention, which produces many spots at arbitrary 
positions deviating from lattice points, requires far smaller resolution and can satisfy the rigorous 
requirement. 
[Fig. 4] 

Fig. 4 is an intensity distribution of a CO2 laser beam adopted by Embodiments 1 to 4 in 
a two-dimensional representation. The intensity distribution is assumed to be a Gaussian 
distribution. 
[Fig. 5] 

FIG. 5 is a DOE pattern designed by the teaching of Embodiment 1 of the present 
invention for producing diffraction spots (Pi= 3.93701 mrad) of Fig. 3. Thickness variation is 
represented by white-black tone variation. There is no repetition of partial patterns in the DOE. 
There is no concept of unit patterns which are identical small pattems. 
[Fig. 6] 

Fig. 6 is a picture of diffracted beam spots on an image plane, which are produced by 
diffracting a CO2 laser beam by the DOE designed by Embodiment 1, wherein a side of the 
diffracted square pattern is 39.75 mrad. Spots appear not on exact lattice points but on distorted 
lattice points which are similar to the deformed lattice pattern of Fig. 3. The beam intensity of 
the spots is uniform and the power of the divided beams is uniform. 
[Fig. 7] 

Fig. 7 is another DOE pattern designed by Embodiment 2 for producing another object 

diffracted beam pattern of the Fig. 3 pattern (Pi = 39.37008 mrad) multiplied by 10. Since the 

39 



object pattern is ten times as large as the Fig. 3 pattern, the size of the pattern of the DOE is 
reduced to one tenth of the first DOE of Fig. 5. Like Fig. 5, the thickness distribution of the 
DOE is represented by variation of white-black monochromatic tones. 
[Fig. 8] 

Fig. 8 is a picture of diffracted beam spots on an image planet, which are produced by 
another DOE designed by Embodiment 2 for making an object diffraction pattern (Pi= 39.37008 
mrad). The size of the diffracted beam picture is 331.25 mrad. Fluctuation of beam spot 
positions seems to be smaller than the previous diffraction pattern of Fig. 6 (Embodiment 1). 
Fig. 8 shows that the DOE satisfies the purpose of reviving the object diffraction pattern. The 
reason why the fluctuation is reduced in Fig. 8 is that the ratio of 5i and 82 to P, and P2 is one 
tenths (1/10) of Embodiment 1. The size of the spots is diminished because the ratio of 6, and 62 
to Pi and P2 are reduced to 1/10 of Embodiment 1 similarly. 
[Fig. 9] 

Fig. 9 is a pattern of a DOE of Embodiment 3 which is made by unifying a Fraunhofer 
type DOE with a converging lens for producing an image beam pattern with fluctuation (position 
deviation) like Embodiment 1 . Thickness distribution of the converging DOE of Embodiment 3, 
which is Embodiment 1 (Pl=3.93701 mrad) plus the converging lens, is represented by white- 
black varying monochromatic tones. There is no square unit pattern in the DOE. Concentric 
patterns are obscurely hidden in the DOE. Superposition of the converging lens produces 
implicit concentric fringes. Pitches of repetition of the DOE are large, since object spots have 
narrow intervals. 
[Fig. 10] 

Fig. 10 is a pattern of another Fresnel type DOE of Embodiment 4 which is also made 

by unifying a Fraunhofer type DOE with a converging lens for making a beam spot group which 

deviates fi-om lattice points like Embodiment 2. Namely, the DOE of Embodiment 4 is built 
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with the Embodiment 2 DOE (Pi= 39.37008 mrad) and the converging lens. Thickness 
distribution of the Embodiment 4 DOE is represented by variation of white-black 
monochromatic tones. There is no unit pattern in the DOE. Superposition of the lens vanishes 
unit patterns and makes the concentric patterns. Pitches of repetition of the DOE pattern are 
small, since object spots have wide intervals (Pi = 39.37008 mrad). 
[Description of Symbols] 
[0164] 

C: pixel (picture element) 
T : unit pattern 

DOE: diffractive optical element 

FFT: the Fast Fourier Transformation 

a: x-direction size (width) of a pixel 

b: y-direction size (length) of a pixel 

P: number of pixels aligning in an x-direction 

Q: number of pixels aligning in a y-direction 

p: x-direction diffraction angle order of beams diffracted by a DOE just on regularly aligning 
lattice points (discrete) 

q: y-direction diffraction angle order of beams diffracted by a DOE just on regularly aligning 
lattice points (discrete) 

Wpq: complex amplitude of a beam diffracted in a (p, q) direction by a DOE (zero for other than 
integers (p, q)) 

Ipq: intensity of a beam diffracted in a (p, q) direction by a DOE 

k: number of diffracted beams ' 

mic: diffraction order of the k-th diffracted beam 

tmn: amplitude transmittance of the (m, n)-th pixel of a DOE 
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A,: wavelength of a laser beam 
A: size of a unit pattern on a DOE (= aP, bQ) 
P : x-direction pixel number in a unit pattern 
Q: y-direction pixel number in a vuiit pattern 

a: x-direction diffraction angle of a beam diffracted by a DOE (continual value) 

p: y-direction diffraction angle of a beam diffracted by a DOE (continual value) 

W(a, p): complex amplitude of a beam diffracted in a (a, p) direction by a DOE (continual value) 

I(a, p): intensity of a beam diffracted in a (a, p) direction by a DOE (continual value) 

U: x-direction angular resolution of diffracted beams 

V: y-direction angular resolution of diffracted beams 

R: x-direction aligning pixel number of a DOE 

S: y-direction aligning pixel number of a DOE 

G: x-direction aligning unit pattern number in a DOE (=R/P) 

H: y-direction aligning unit pattern number in a DOE(=S/Q) 

E: merit function 

Wi. W2, W3: weights multiplied to squares of deviations of parameters in a merit function 

R': x-direction aligning pixel number of a DOE which is enlarged from an original DOE by 

adding vacant pixels (t„„=0) at peripheries 

S': y-direction aligning pixel number of a DOE which is enlarged from an original DOE by 
adding vacant pixels (t,„„=0) at peripheries 



42 



